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Entanglement is the quintessential quantum mechanical phenomenon understood to lie at the 
heart of future quantum technologies and the subject of fundamental scientific investigations. Mix- 
ture, resulting from noise, is often an unwanted result of interaction with an environment, but is 
also of fundamental interest, and is proposed to play a role in some biological processes. Here we re- 
port an integrated waveguide device that can generate and completely characterize pure two-photon 
states with any amount of entanglement and arbitrary single-photon states with any amount of 
mixture. The device consists of a reconfigurable integrated quantum photonic circuit with eight 
voltage controlled phase shifters. We demonstrate that for thousands of randomly chosen configura- 
tions the device performs with high fidelity. We generate maximally and non-maximally entangled 
states, violate a Bell- type inequality with a continuum of partially entangled states, and demonstrate 
generation of arbitrary one-qubit mixed states. 



Quantum mechanics is known to allow fundamentally 
new modes of information processing [TJ |2] , simulation 
[3j 2], and communication [5 as well as enhanced pre- 
cision of measurement and sensing [6] [7]. Single pho- 
tons provide a particularly promising physical system 
with which to develop such quantum technologies [8]— 
due to their low noise, high speed transmission and ease 
of manipulation at the single photon level — and have long 
been a leading approach to exploring fundamental quan- 
tum science. The ability to precisely prepare, control and 
measure multi-photon states therefore holds considerable 
scientific and technological interest. 

Recently it has been shown that it is possible to 
miniaturize quantum optical circuits using optical fibre 
[9] [10] and integrated waveguide chips [TTHT7I . Mono- 
lithic waveguide circuits are inherently stable and can be 
many orders of magnitude smaller than their bulk optical 
equivalents, enabling the fabrication of multi-purpose, re- 
configurable quantum circuits of unprecedented size and 
complexity. Control of a single phase shifter in this archi- 
tecture has been used to manipulate time bin qubits [18] 
and up to four photons in two spatial modes [12], how- 
ever, the large-scale reconfigurability required to gener- 
ate arbitrary multi-photon states, including mixture and 
entanglement, has so far been out of reach. 

Here we report an integrated quantum photonic de- 
vice comprised of a two-qubit entangling gate, several 
Hadamard-like gates, and eight variable phase shifters, 
and demonstrate that it can be reconfigured with high 
fidelity across the complete space of possible configura- 
tions. We use this device to generate all four Bell states 
and perform quantum state tomography on them, to re- 
alise a Bell inequality "manifold" — obtained from a con- 
tinuum of measurement settings and states with a vari- 
able amount of entanglement — and to prepare and mea- 
sure arbitrary single-photon mixed states. 

A reconfigurable quantum photonic circuit 

The device described here is a silica-on-silicon entangling 
circuit, shown in Fig. [I] Two photonic qubits A and B 
are encoded in pairs of waveguides — path or dual rail 
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FIG. 1: A two-photon reconfigurable quantum circuit for gen- 
erating, manipulating and detecting entanglement and mix- 
ture, (a) Quantum circuit diagram consisting of pairs of 
Hadamard-like gates H' = e i7T/2 e- ilT(7z /4 H e - ilT(7z /4 (where H 



is the usual Hadamard gate) and Rzi^) — e 
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that together implement Uij(4>j, 4>k), and two H' gates and a 
controlled-SiGN or cz gate, that together implement a CNOT. 
(b) Waveguide implementation of the circuit composed of di- 
rectional couplers and voltage controlled thermo-optic phase 
shifters drawn as A. Directional couplers with splitting ratio 
r\ — 1/3 are marked •, all other couplers have r\ — 1/2. 



encoding. These two qubits are input in the logical zero 
state |0a) ® |0#) — i.e. a single photon in each upper 
waveguide — and are then acted upon by the quantum 
circuit shown in Fig. [I] 

The first part of this circuit enables arbitrary state 
preparation of each qubit. The central part of the 
circuit implements a maximally entangling postselected 
controlled-NOT (cnot) logic gate pQ — the canonical two- 
qubit entangling gate. The cnot gate is a postse- 
lected linear optical gate that works with probability 1/9 
[T9] [20] . The final stage of the circuit is the mirror image 
of the first stage and is followed by measurement in the 
computational basis, which together enables projective 
measurement of each qubit in an arbitrary basis. 

The initial and final stages of the device can be re- 
configured, and are each implemented using two MZ in- 
terferometers, each composed of two voltage-controlled 
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FIG. 2: Classical and quantum interference fringes, (a) Interference fringe measured at the two outputs of a single Mach-Zehnder 
(MZ) interferometer on the chip. Experimental data are presented as black circles. Solid lines show fits, (b) Hong-Ou-Mandel 
dip, measured using a single MZ interferometer as a beamsplitter. Two-photon coincidence counts are shown as black circles. 
The red line shows a fit to this data with Gaussian and sine components, due to quantum interference and determined by the 
spectral filters used, and a linear term accounting for slight decoupling of the source. The blue line shows a fit to the measured 
rate of accidental coincidences, with Gaussian and linear components. Error bars in both figures assume Poissonian statistics. 



thermal phase shifters and two directional couplers [T2] . 
This architecture allows reconfigurable single-qubit uni- 
tary operations to be performed: In general any uni- 
tary in SU (2) can be realised using three phase shifters 
and an MZ interferometer [21] as Uarbi^a^b^c) — 
e iif c a z /2 e iif b a y /2 e iif a a z /2 _ jj ere we use two phase shifters 

per MZ to realise Ui(ipb,<p c ) = e -^ cCT */ 2 e"^ b ^/ 2 at the 
input, and Uf = llj at the output, of each qubit. This 
is adequate for arbitrary state preparation and measure- 
ment, up to a global phase. Explicitly, the entire circuit 
shown in Fig. [I] implements the unitary matrix 



£7/(05,06) ®U S (07,0s) 



' UcNOT ' 



Ui (01,02) ® Ui (03,04) 



where U C not = |00> (00| + |01) <01| + |11) <10| + |10) (11| 
and 0i_8 are set by the external control voltages. 

Benchmarking of reconfigurability 

This circuit can be reconfigured to perform a number of 
different tasks, including arbitrary two-qubit pure (en- 
tangled) state preparation, arbitrary one-qubit (mixed) 
state preparation, state tomography, process tomogra- 
phy, etc., as detailed below. 

In order to characterize the precision and accuracy 
with which the device can be reconfigured we injected sin- 
gle photons into the device via a polarization maintain- 
ing optical fibre array, and measured interference fringes 
across each of the eight phase shifters on the chip, finding 
an average contrast C = 0.988 ± 0.008. (See Appendix 
for details.) An example of a pair of such fringes is shown 
in Fig. [2^i . From these measurements, we estimate the 
average accuracy in phase across all eight heaters to be 
5^ rsj 0.05 radians. 

In addition to high-fidelity classical interference, as 
demonstrated in Fig. [2^i, the cnot gate in the mid- 
dle of the circuit shown in Fig [I] relies on high-fidelity 



quantum interference [16 j. Fig. [2]o shows a Hong-Ou- 
Mandel dip [22] measured at a single MZ interferometer 
(that containing 0i) on the chip. We produced degen- 
erate photon pairs, sharing the same spectral and polar- 
ization mode, via type-I spontaneous parametric down- 
conversion [23] (see Appendix) which were injected into 
the chip as shown in Fig. [I] The phase in the interfer- 
ometer was then set to 7r/2, rendering it equivalent to a 
1/2 reflectivity beamsplitter, and the two-photon coinci- 
dence count N across the outputs of the interferometer 
was measured as a function of an off-chip optical delay 
between the arrival times of the two photons. The vis- 



ibility of the dip V = (N a 



quantum 



)/N c , 



was measured to be 0.978 ± 0.007, taking into account 
the measured rate of accidental coincidences [24]. Bire- 
fringence, mode mismatch, and other such imperfections 
in the circuit could limit the visibility of this dip. The 
high visibility of the HOM dip therefore indicates the 
high quality of the device. 

Having observed high-fidelity classical and quantum 
interference at individual MZ interferometers on the chip, 
we then used a stochastic method to characterise the op- 
erational performance of the quantum circuit as a whole, 
across the full space of possible configurations. We chose, 
at random, 995 vectors ipj representing possible configu- 
rations of the device 



<Pj = 0i>02>--->08 



(1) 



with < 0^ < 27r. Injecting photon pairs as before , the 
probability-theoretic fidelity / = ^ k ^/pk • p' k between 
experimentally measured coincidence probabilities at the 
output of the device (poo, Poi, Pio, Pn) an d the ideal 
theoretical values (p' 00 , £>o l5 Pioi Pii) was calculated for 
each (fj. The statistical distribution of these fidelities 
is shown in Fig. [3^i. The average fidelity across 995 
configurations (equivalent to many truth tables in many 
bases) was measured to be 0.990±0.009 with 96% of con- 
figurations (fj producing photon statistics with / > 0.97. 
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FIG. 3: Statistical fidelity of photon coincidence count rates. 
The histogram shows the distribution of statistical fidelity be- 
tween ideal and measured coincidence count rates, over 995 
sets of eight randomly selected phases (p. 96% of phase set- 
tings produced statistics corresponding with theory to / > 
0.97. 





FIG. 4: Real parts of the density operators of the Bell states 
|$ + ), |$~), |^ + ) and |^~) ( (a) , (b) , (c) , (d) respectively), 
generated and characterized on-chip. 



This result depends on simultaneous high fidelity quan- 
tum and classical interference, as well as accurate and 
precise joint control of all eight phase controllers. Poor 
performance of any of these component parts would re- 
sult in lower fidelity output for some subset of {ipj}. The 
high fidelity operation observed in these tests of reconfig- 
urability bodes well for the operations described below. 

Generating and characterising entanglement 

Entangled states of quantum systems are the fundamen- 
tal resource in quantum information and represent the 
most nonclassical implication of the formalism of quan- 
tum mechanics. The circuit shown in Fig. [I] can be used 
to prepare a continuum of entangled, partially entangled, 
and separable states with only computational-basis prod- 
uct states as input. 

In order to demonstrate this ability, we first prepared 
and analysed each of the four maximally entangled Bell 
states. Inputting the \0a) \0b) state as before, the state 
preparation stage of the circuit was used to generate each 
of the superposition states I 0b), \^a) | 1b), where 
|±) = (|0) ± |1)) /y/2 , at the input of the CNOT gate. 
The corresponding Bell states (| < ^ )± ) and |^ r± ) respec- 
tively) are ideally produced by the CNOT gate. 

We used the arbitrary single-qubit measurement capa- 
bility of the circuit to perform maximum-likelihood quan- 
tum state tomography (QST) [25 on these four states: 
phase shifters 05-8 were used to implement each of the 16 
measurements necessary to reconstruct the density oper- 
ator of the state. The measured density matrices of all 
four Bell states are shown in Fig. |4j with quantum state 
fidelities F = (Try/^p^p exp ^p^) 2 of 0.947 ± 0.002, 
0.945 ± 0.002, 0.933 ± 0.002, and 0.885 ± 0.002 respec- 
tively. 

A Bell-type inequality manifold 

The Clauser, Home, Shimony, and Holt (CHSH) [26] test 



of local hidden-variable models of quantum mechanics 
requires that the sum 

S = (ixBi) + {A X B 2 ) + {A 2 B 1 ) - (A 2 B 2 ) (2) 

satisfies the Bell-CHSH inequality — 2 < S < 2 for any 
local hidden- variable model, where B{ are measure- 
ment operators chosen by two observers, Alice and Bob. 

The Bell-CHSH experiment provides a well-known test 
for the presence of entanglement that we use here to ex- 
amine the performance of the device, as it is reconfigured 
across a large parameter space. Specifically, we use 4>\-a 
and the CNOT gate to prepare the state 

IVW> = ±= [(1 - e"*) |00) + (1 + e te ) |11)] , (3) 

where a = (pi. By changing a it is thus possible to tune 
continuously between two orthogonal maximally entan- 
gled states: for a = 0,7r, |Vw) is a product state, and 
with a = 7r/2,37r/2, |?/w) is the maximally entangled 
state ^= (| 00) ± i |11)) (up to a global phase). In the 
course of this preparation, we pass through a continuum 
of partially entangled states. In order to evaluate S we 
make four two-qubit measurements on the state emerging 
from the CNOT gate, which correspond to combinations of 
observables chosen by Alice and Bob. While Alice's two 
measurement settings, <j)Q = 7r/4, — 7r/4, do not change, 
Bob's two measurement settings are varied continuously, 
as 08 = P, ft + tt/2. We measured S(a, (3) for a G [0, 2tt] 
and j3 G [0, 27r], with step size 27r/15, producing the "Bell 
manifold" shown in Fig. [5j We measured maximum and 
minimum values of S of 2.49 ± 0.03 and -2.54 ± 0.03 
respectively. Errors were again determined by a Monte- 
Carlo technique, assuming Poissonian statistics. The co- 
efficient of determination 

p2 _ -| T i) 2 ( a\ 

R - 1_ (4) 
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FIG. 5: CHSH manifold, (a) The Bell-CHSH sum S, plotted as a function of phases a and f3. In the a axis, the state 
shared between Alice and Bob is tuned continuously between product states at a = 0,7r and maximally entangled states at 
a = 7r/2, 37t/2. The f3 axis shows S as a function of Bob's variable measurements, which can be thought of as two operator-axes 
in the real plane of the Bloch sphere, fixed with respect to each other at an angle of 7r/2 but otherwise free to rotate with angle 
P between and 2tt. The blue curves show a projection of the manifold onto each axis. Yellow contours mark the edges of 
regions of the manifold which violate — 2 < S < 2. Red lines on the axes also show this limit, (b) Experimentally measured 
manifold. Data points are drawn as black circles. Data points which violate the CHSH inequality are drawn as yellow circles. 
The surface shows a fit to the experimental data. 



where Si are experimentally measured values of the Bell- 
CHSH sum, S is the average over Si, and Ti are the 
ideal theoretical values of 5, evaluated across the whole 
manifold, is 0.935. Here the ideal value is one, the worst 
case yields zero. 

Generating and characterising mixture 

Mixture is often associated with noise or decoherence in 
quantum processes and its deliberate and controlled im- 
plementation is critical for characterisation of devices; 
furthermore it has been shown that quantum comput- 
ing can be performed despite mixture [27] . More signifi- 
cantly, recent work has suggested that decoherence may 
play an important role in biological processes that ex- 
hibit quantum coherence [28j [29]; photonic waveguide 
systems show great promise for simulating these pro- 
cesses [30j[3T], however, such simulations will require the 
controlled introduction of mixture. 

By tracing over one of the two output photons, our 
device can prepare an arbitrary state of a single qubit 
including any amount of mixture. The amount of entan- 
glement in the two qubit state prepared by the first stage 
of the device determines the degree of mixture, which can 
range from zero (a pure state) to one (a maximally mixed 
state) . 

In general the state 

Wont = <*7 \Qa0b) +a5 \0 A 1 B ) + /3 7 \1 A 1 B )+P6 \1 a b ) 

(5) 

is generated after the cnot gate in the circuit, where a, 
/3, 7, 5 are complex parameters related to 0i_4. Trac- 
ing out the second qubit, we find the reduced density 



operator of qubit A, 

PA = M 2 |o><o| + a/r(7<r + <h*)|o><i| (6) 

+ /?a*( 7 r + 5 7 *)|l)(0| + |/?| 2 |l)(l|. 

By choosing a, /3, 7, via setting 0i_4, the amount of 
mixture in this reduced density matrix can be continu- 
ously varied between and 1. 

We chose 119 target states with various amounts of 
mixture, at random by the Hilbert-Schmidt measure [32] , 
then generated each state and reconstructed its density 
matrix by maximum likelihood state tomography using 
phase shifters and 0g. Fig. [6]shows the fidelity of these 
reconstructed states. The average quantum state fidelity 
across all 119 states was measured to be 0.98 ±0.02, with 
91% of states having fidelity > 0.95. We then chose 63 
specific mixed states that mapped out the symbol 
inside the Bloch sphere, and generated them with high 
fidelity (Fig. [6| inset). 

We note that the device shown in Fig. 1 could be 
used to generate mixture by applying random voltages 
to phase shifters on a single qubit, without the need for 
entanglement. We chose to use the entanglement ap- 
proach as a more demanding test of our device, demon- 
strating sufficient control to obtain the data shown in Fig. 
6. An advantage of using the entanglement approach in 
practical applications is that it does not require pseudo- 
/quantum-random number generators. 

Discussion 

Quantum information science and technology with pho- 
tons will require circuits that are complex, stable and 
highly reconfigurable in a straightforward manner. High 
fidelity production and measurement of states of arbi- 
trary entanglement and mixture will be essential for char- 
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FIG. 6: Histogram showing the statistical distribution of 
quantum state fidelity between 119 randomly chosen single- 
qubit target states and the corresponding mixed states gen- 
erated and characterized on-chip. Inset: \I/ drawn in the the 
Bloch sphere using 63 mixed states, again generated and char- 
acterized on-chip. These states are chosen from the real plane 
of the sphere for clarity. Note that each point is derived from 
a different bipartite partially entangled state. 

acterisation of quantum devices, and will provide a re- 
liable means to test the unique properties of quantum 
physics. The generation of mixed states may also be im- 
portant in quantum photonic analogues of biochemical 
systems that rely on decoherence [28] [29] . 

Although on-chip polarization encoding is possible 
[15], the inherent interferometric stability of integrated 
optics makes path encoding of qubits a natural choice, 
with the further advantage that encoding of higher- 
dimensional qudits [33] is immediately possible. This 
is in contrast with bulk optics, where two- level polar- 
ization encoding is more natural, and stable path en- 
coding requires a considerable resource overhead. Fur- 
thermore, this architecture could be used to manipulate 
hyper-entanglement [34] encoded with multiple degrees of 
freedom [35, 36 . Circuits such as the one presented here 
could be used in conjunction with adaptive (classical) al- 
gorithms to bypass the need for calibration of the phase 
shifters in particular applications. For example, repeated 
measurement and feedback onto the voltage-controlled 
phase shifters based on comparison of the output state 
with a desired target state could be used to reconfigure 
the circuit via a genetic algorithm. 

Appendix 

Device: The waveguide device was fabricated on a sil- 
icon wafer, upon which a 16/im layer of undoped silica 
was deposited to form the lower cladding of the waveg- 
uides. 3.5-/im-wide waveguides were then patterned in a 
3.5-/im layer of silica doped with germanium and boron 
oxides. A 16-/im layer of silica, doped with phospho- 
rous and boron so as to be index-matched with the lower 
layer, constitutes the upper cladding. Resistive heaters 
and corresponding electric contacts were then patterned 



in metal on top of the chip using standard lithographic 
techniques. Dimensions of the chip are 70mm x 3mm. 

Photon Source: Degenerate pairs of 808 nm photons were 
generated by focusing a 404 nm, 60mW laser onto a 2mm 
thick, Bisumuth Borate BiB406 (BiBO) nonlinear crys- 
tal, phase matched for type I spontaneous parametric 
down conversion (opening angle of 3 degrees). 

Photon pairs were spectrally filtered using 3nm full- 
width at half maximum interference filters. The in- 
terference niters were designed with central wavelength 
808 nm and were tilted to ensure that photons pairs 
were identically filtered. Photons were then collected 
into polarisation maintaining fibre (PMF) using 11mm 
aspheric lenses. Typical two-photon coincidence count 
rates of 100kHz were achieved using ~ 60% efficient, sili- 
con based avalanche photo-diode single photon counting 
modules (SPCM). The photons were then launched into 
the waveguide chip by butt-coupling arrays of PMF with 
250/im spacing (matched to the input and output waveg- 
uide pitch). Photons were collected from the output of 
the chip also using arrays of polarisation maintaining fi- 
bre and detected with fibre coupled SPCMs. A typical 
facet-to-facet coupling efficiency of ~ 60% was achieved. 

Calibration of Phase Shifters: Each thermal phase shifter 
4>i has a nonlinear phase voltage relationship: 

MK) = en + fcV 2 + 7 ^ 3 + S % V\ (7) 

where V{ is the voltage applied across phase shifter i and 
(j>i is the resulting phase shift, ctj, ft, 7^ and Si are 
real numbers associated with the response of a partic- 
ular heater. Each phase shifter can be seen to occupy 
one particular MZ interferometer in the circuit. 02 and 
05 can be seen as acting on a single, lossy MZ interfer- 
ometer. 

Each phase shifter in the circuit was calibrated as fol- 
lows: Bright light from an 810nm 1.3mW laser was in- 
jected into one input port of each MZ, and the intensity 
at each output port was measured as a function of the 
voltage applied across the heater, which was swept lin- 
early between 0V and 7V. This produced classical inter- 
ference fringes, distorted by the nonlinear phase-voltage 
relationship ([7]). We then fitted the function 

I(V)=A(l-C-cos 2 (4>(V)/2)) (8) 

to each set of experimental data with A, C, a, ft 7 and S 
as fitting parameters. This yields the complete (approxi- 
mate) phase-voltage relationship for each heater. No evi- 
dence of crosstalk between phase shifters (due to thermal 
effects or otherwise) was observed in these experiments. 

Single-Photon Fringes: We measured single photon 
fringes for each phase shifter using photon pairs from the 
source. Injecting photons from one arm of the source into 
the chip, we counted coincidences between single photon 
events from a particular output of the interferometer, and 
those from the other arm of the source. This approach 
largely mitigates the contribution of SPCM dark counts. 
The fit shown in Fig. [2]has the form A(l - C • cos 2 (0/2)) 
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where C is the fringe contrast and A is the amplitude. 
Each fringe is normalized with respect to its amplitude 
for clarity. 
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